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Abstract. We studied the ISM distribution in and around the star forming cloud L1251 with optical star counts. 
A careful calculation with a maximum likelihood based statistical approach resulted in B, V, R, I extinction 
distributions from the star count maps. A distance of 330 ± 30pc was derived. The extinction maps revealed an 
elongated dense cloud with a bow shock at its eastern side. We estimated a Mach number of M ~ 2 for the bow 
shock. A variation of the apparent dust properties is detected, i.e. the Rv ~ Av / Eb-v total to selective extinction 
ratio varies from 3 to 5.5, peaking at the densest part of L1251. The spatial structure of the head of L1251 is well 
modelled with a Schuster-sphere (i.e. n=5 polytropic sphere). The observed radial distribution of mass fits the 
model with high accuracy out to 2.5pc distance from the assumed center. Unexpectedly, the distribution of NH3 
1.3 cm line widths is also well matched by the Schuster solution even in the tail of the cloud. Since the elongated 
head-tail structure of L1251 is far from the spherical symmetry the good fit of the linewidths in the tail makes it 
reasonable to assume that the present cloud structure has been formed by isothermal contraction. 
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1. Introduction 



The Cepheus Flare is an extended complex of molecular 
clouds, where luminous stars are clearly absent, and it 
harbors dark clouds forming low mass s tars (for details 
see the comprehensive study o f lKunI lll99?^ on this region) . 
L1251 (^2000 = 22''36™l;(52ooo = +75°16') is a dark cloud 
in th is region at a distance of 300 ± 50pc ijKun fc Prustil 
I1993 I) and apparently belongs to this complex. Already at 
the advent of molecular radio astronomical studies in the 
late sixties the cloud was detected am ong the ten brightest 
OH em ission sources on the sky by 'Cud aback fc Heilel 
l)l969(l . The cloud was also listed aniong sources of stro ng 
formaldehyde emission llDieteilll973l:ISume et al.lll97,'il) . 

After these early successes of molecular radio obser- 
vations LI 251 apparentl y escaped t he attention of ra- 
dio observers. Although iKunI lll982t) discovered a num- 
ber of Ha emission stars associated with L1251, only ap- 
peared again in mole c ular radio studies i n the late eighties 
jSato fc FukTil Il989t IZhou et alJ Il989l iBenson fc Mver j 
Il989|) . IzTiQiTetalJ l)l98 i!l) recognized an ammonia core in 
the dense part of the cloud and comparing the large num- 
ber of Ha objects associated with the cloud . ISato fc Fuku] 
l)l98 9^ concluded that the star formation efficiency is 
anomalously high in L1251. From all of these facts it 



was obvious that the star formation processes in this re- 
gion needed further detailed investigations. This moti- 
vated Kim & Prusti (1993) to study the properties and 
distribution of faint IRAS point sources along with the Ha 
emission objects. Their study indicated that L1251 has 
been forming low-mass stars with an efhciency higher than 
usually encountered in dark clouds. The eastern, head re- 
gion of the cloud has been found to contain more evolved 
YSOs than the western (tail) side. 

A census o f den se cores carried out by 
iToth fc WalmslevI l|l996|) based on the NH^ 1.3 cm 
line using the 100 m dish at Effelsberg resulted in 
detection of eight ammonia cores with typical size of 
FWHM « 2' (0.2pc at 350pc distance). Five of the 
cores were found to be gravitationally bound. L1251 
was recently surveyed in several mm lines, the struc- 
ture and kinematic s of the co r es we r e studied e .g. hy 



ICasdh et ahl lEonj. IPee et a,1.l IPee et a,1.l I^M 

and lNikolic et al.l l|2003() . Both infall and outward motions 



we re detected i n the gravitationally unstable cores. 

iToth et "all l)l995|) tried to explain t he shape of the 
cloud, as seen in the CO observations ijGoodman et alJ 
Il993t ISato et al]ll994|) , by hydrodynamical modelling as- 
suming an encounter of a dense molecular cloud with an 
external shock, probably originating from a nearby super- 
nova explosion. The reality of this assumption is supported 
by the detection of a soft X-ray excess region east of the 
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cloud byl&i nier et aD l)l989() . Further confirmation of an 
SN explosion co mes from the spa, ce motion of the runaway 
star HD203854 llKun et al.ll2f)n(il) . 

Although the dense gas component of the cloud is 
well studied, much less is known about the distribution 
and prop erties of inters tellar dust in and around LI 251. 
Recently iKandori et al.l l(2003il studied the extinction of 
the dust component using B, V, R, /star counts and found 
high Rv values indicating grain growth in the head of the 
cloud. 

In the present paper we study the spatial distribution 
of dust, the mass, and the basic physical properties in the 
cloud by means of new optical extinction maps in B, V, 
R, /colors. In a subsequent paper the thermal emission of 
the dust grains in the dense regions of the cloud will be an- 
alyzed on the basis of far- infrared (120 and 200/LtTO) maps 
obtained from the Infrared Space Observatory (Toth et 
al. in prep.). The two complementary studies will clarify 
the spatial distribution of dust with respect to the distri- 
bution of the dense gas, and indicate how the star forma- 
tion activity and the encounter with the external shock 
has modified the properties of the dust particles. To ob- 
tain direct information on the spatial distribution of dusty 
material the study of star counts is still one of the most 
reliable approaches. The basic aim of this paper is to carry 
out such analysis. 

2. Input data 

To study the surface distribution of the optical extinction 
we obtained star counts in B, V, R, I colors based on pho- 
tographic observations with the 60/90/180 cm Schmidt 
telescope of the Konkoly Observatory. The plates were 
digitized with a pixel size of 20 microns in a 1.5 x 1.5 
degree field around LI 251 using the PDS microdcnsito- 
meter of the Vienna Observatory in 1991. We scanned 4, 

3, 3, and 2 plates in B, V, Rj, Ij colors, respectively. The 
scans were processed with the ROMAPHOT photometric 
programme integrated into the MIDAS data analysis pack- 
age. The plates were calibrated via CCD observations per- 
formed with the 1.23 m Ritchey-Chretien tel escope of the 
Ger man-Spanish Observatory, Calar Alto 1 Balazs et alJ 
Il992^) . The limiting magnitude of the photographic sur- 
vey was 19.0, 18.5, 17.5, and 16.5 in B,V,Rj, and Ij, 
respectively. Although these figures are less than could 
be obtained from the Digitized Sky Sur vey maps (see e.g. 
the extinction maps of ICambresvl 1 1999^ ) . the well defined 
color system, the much better photometric calibration and 
four colors, however, support using our data. The esti- 
mated completeness of star counts was about 1.5 magni- 
tude above the detection limit of the plates. The final star 
count maps of the region were obtained by counting the 
stars in each color on pixels of 6' x 6' size and a 2' mesh of 
the star count maps was selected in both directions. (This 
resolution approxi mately corr espon ds to those of the 100 
lim IRAS maps). iKunI l)l982|) and iKun fc Prustil l|l993|) 
lists 12 Ha objects (candidate pre-main sequence stars) 
apparently associated with the cloud. Their effect on the 



R and / star counts might be significant. Consequently, 
we have omitted them from further analysis. 

3. Extinction maps from star counts 

Several studies indicated empirically that the surface dis- 
tribution of st ar counts i s an e x cellent tracer of optical e x- 
tinction. fe.g. iDickinaiil l)l978|) : ICernicharo et all l)l985|) ). 
This means that there is a simple linear relationship be- 
tween the logarithmic star counts and the extinction 
of a dust cloud: 



ad ^ a X log{N{m)) + b{m) 



(1) 



where N{m) is the cumulative star count up to a given m 
limiting magnitude, the a constant and h{m) depend on 
the Galactic longitude and latitude; —h/a = log(A'o(m)) 
measur es the logarithm ic star count in an extinction-free 
region |Dickmanlll978l) . 

In the following we calibrate this expression using a 
multivariate statistical method, the fc-means clustering, 
and a maximum likelihood procedure. Based on this cal- 
ibrated relationship we assign extinction values to each 
star count pixel yielding an extinction map of the cloud. 
The procedure gives as a byproduct the distance of L1251 
which enables us to calculate the mass of the cloud. Since 
we derive the extinction of the cloud in different colors we 
also discuss the ratio of the selective to total extinction. 

3.1. Definition of areas of similar extinction 

The B, V, R, /star count data on the pixels defined above 
represent a distribution of points in a four dimensional 
B, V, R, I parameter space. Following the relationships 
between star counts and extinction we assumed that the 
regions of equal extinction have equal star counts on the 
maps. In other words, looking for areas of equal extinction 
means searching for points of similar coordinates, i.e. ly- 
ing close to each other in the four dimensional parameter 
space, made up of the pixel values of the star count maps. 

The linear relationship between the logarithmic star 
counts and interstellar extinction predicts a one dimen- 
sional manifold in the four dimensional parameter space, 
stretched by the extinction, and there is a Poisson noise 
superimposed on it by the star counts. Therefore, to con- 
vert star counts into optical extinction we had to divide 
the whole star count map into regions of equal counts. 
According to Equation the logarithms of the star 
counts scale linearly with the extinction; therefore we used 
the logarithmic B, V, R, I star counts in this procedure. 

In order to divide the points into groups of similar 
extinction in the parameter space we invoked the tech- 
nique of fe-means clustering (see e.g. iMurtagh fc HeclJ 
l)l987l) ). K-means clustering orders the points in the pa- 
rameter space into k groups. The k number of the groups 
should be specified before running the clustering proce- 
dure. Assigning the points to any of the groups proceeds 
on the basis of some distance measure between the points. 
We used squared Euclidean distance. 
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L1251 : surface distribution of clusters 
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Fig. 1. Result of the clustering, reprojected onto the ce- 
lestial sphere. The areas of same grey level represent pixels 
in the map belonging to the same cluster and have similar 
extinction. We displayed the cloud in rectangular coordi- 
nates in a plane tangential to the celestial sphere. North 
is at the top and East is to the left. The (0,0) position 
corresponds to azooo ^ 22''36™1; (^zooo = +75°16'. 

Table 1. Summary of the star count analysis of L1251, 
number of stars in each cluster and in each color 



FIELD 


AREA{SqD) 


Nb 


Nv 


Nr 


Ni 


1. 


0.9880 


3982 


3929 


3939 


2723 


2. 


0.8770 


1883 


2045 


2186 


1530 


3. 


0.1515 


101 


143 


146 


138 


4. 


0.0605 


7 


13 


14 


15 


5. 


0.0090 


1 


2 


4 


4 



There are no definite criteria for fixing the value of 
k. By trial we selected k = 5 enabling clear separation 
of the high extinction regions from those of low extinc- 
tion and ensuring enough stars in each cluster for reliable 
analysis. Figure^shows the result of the clustering, repro- 
jected onto the celestial sphere. Areas of same grey level 
in the map represent pixels belonging to the same clus- 
ter and consequently, having similar extinction. Table ^ 
summarizes the number of stars in each subregion (clus- 
ter) in each color. One can infer from the data of this 
table that the 5th cluster is scarcely populated, therefore 
we excluded it from further analysis. 

The central, densest part of L1251 is clearly separated 
from the surrounding lower extinction region. It resem- 
bles a bullet moving at supersonic speed with respect to 
the ambient medium. The less dense area surrounding the 
bullet-like main body of the cloud has a form of a bow 



shock. Accepting this view one can estimate the speed of 
the cloud relative to the ambient medium. The interaction 
of the cloud with its surroundings is probably the key is- 
sue in understanding the history of star formation. We will 
return to this problem in Sect.^ 

3.2. Modelling the effect of extinction on the star 
counts 

3.2.1. Modelling differential star counts 

To calibrate star counts in terms of optical extinction 
we attempted to model the observed apparent magnitude 
distribution of the stars. The s tarting point of our ap - 
proach was the Galactic model of lWainscoat et all l)l992j) . 
Following the basic ideas of this model we assumed that 
the spatial distribution of the stars in the region inves- 
tigated can be satisfactorily described by superposition 
of exponential disks, corresponding to different types of 
stars, and a spheroidal component. The exponential disks 
were defined by their scale heights and local ste llar den- 
sities near the Sun. We used the data of Wainsc oat et alJ 
for the characteristic values of the exponential disks 
in the model. 

To model the effect of the absorbing cloud on the star 
counts we assumed that besides the obscuring matter asso- 
ciated with L1251 there is no other significant dust cloud 
in the line of sight. This assumption is quite reasonable 
due to the high (+1 5 °) gal actic latitude of L1251. As in 
the I Wainscoat et al' I lll992^ model we assumed that the 
diffuse component of interstellar extinction has the form 



(^diffir) = ao sec(6)(l - exp(-rsin(6)/ft,)) 



(2) 



where oq is a constant depending on the color selected, 
b is the Galactic latitude of the cloud and h is the scale 
height of the obscuring material. We added to the extinc- 
tion described above a further component in the form of 
a step function. 



0, if r < rci 
ad, if r > rci 



(3) 



The distance of the cloud, , and the cloud extinction, 
ad, are constants to be estimated within a cluster ob- 
tained by the k-means clustering in the star count param- 
eter space. 

3.2.2. Maximum likelihood estimation of the distance 
and extinction 

The model described in the previous section allows us to 
get the distance and extinction of the cloud using the max- 
imum likelihood (ML) estimation. According to our model 
assumption the probability density of the apparent mag- 
nitude of stars in our observed sample is given by 



A{m\ad,rd) = 
'E.p/o"^.pM'i>("^-51og(r)-5 



(4) 

a{r\ad,rd))r'^dr 
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In the above formula lu is a suitably chosen normalizing 
constant, Dsp{r) is spatial density, <I> is the luminosity 
function assumed to have a Gaussian f orm and sp runs 
over t he spectral types represented in the lWainscoat et all 
model. We used a Monte Carlo simulation to get 
a distribution of apparent magnitudes corresponding to 
the probability density of A{ni\aci,rci)- We used the MC 
simulated data to calculate the numerical values of A for 
the ML procedure. 

The likelihood function in our case can be written as 

n 

L{aci,rci) ^^\og{A{mi\aci,rci)) (5) 

i=l 

where rrii-s are the observed apparent magnitudes in one 
of the colors in our sample. Maximization of L{aci,rci) 
with respect to ad and r^i yields the ML estimation of 
the extinction and the distance of the cloud. 

Performing the ML estimation within all groups given 
by the fc-means clustering and in all colors separately, we 
obtained the results summarized in Tables [21 and |31 The 
results summarized in Table [3 enable us to convert the 
star counts in different colors into extinction. We return 
to this calibration in the following subsection. 

One may use the calculated distance moduli in TableOl 
to get an estimate for the distance of L1251. We computed 
a weighted mean of the data in the table using weights in- 
versely proportional to a. It resulted in a distance modulus 
of 7.58 ± 0.2 corresponding to 330 ± 30 pc. 

3.2.3. Confidence interval for the parameters estimated 

The ML estimation provides a straightforward way to ob- 
tain the confidence interval for the estimated parameters, 
the extinction and distance. Denoting the value of the pa- 
rameters maximizing the likelihood function with a™°^, 
r™°^ and with a*™^ ' ''cl"'^ their true values we have 
asymptotically, if the sample size goes to infinity, 

2[LiaT^, rTl - L{a%^\ r*r )] = xl, k^2 (6) 

In general, k equals the number of parameters esti- 
mated and Xk ^ x-square variable wit h k degrees of free- 
dom (for the proof of this theorem see iKendall fc Stuart! 
l)l976. 1979 )). The probabihty that the true values of the 
parameters are within a certain region in the k dimen- 
sional parameter space is given by 

Pixl < xL) = 1 - (7) 

The projection of this k dimensional domain which is 
given by the xt — Xok inequality yields the confidence in- 
terval of the individual parameter values estimated by the 
ML procedur e. The Xni.,S pairs are tabulated in the liter- 
ature (see e.g iKendall fc StuartI (Il976. 19791) ). The confi- 
dence intervals corresponding to the la levels are given in 
the (T columns of the tables. In the case of Field 4, due to 
the small numbers of stars the confidence interval is not 
closed towards higher extinctions and lower distances. 



Table 4. Linear least squares fitting of the ad — a x 
log(7V(m)) -t- b{m) extinction - star count relation 



color 


a 


o-a 


b{m) 


o-b 


B 


-4.08 


0.23 


5.78 


0.21 


V 


-4.08 


0.16 


5.67 


0.15 


R 


-4.39 


0.30 


5.86 


0.28 


I 


-4.52 


0.18 


5.31 


0.16 



3.3. Star count - extinction conversion 

3.3.1. Verifying the linear log A^(m) - extinction 
relationship 

The results yielded by the ML analysis enabled us to verify 
the star count - extinction relationship in each color in 
our study. Assuming the functional form of the formula 
given by Eq. one can get its constants by a linear least 
squares fitting of the star counts versus extinction given 
in Table El We summarized the results of the least squares 
fitting in Tabled 

Table d clearly demonstrates that the linear relation- 
ship fits the data obtained from the ML analysis. This 
means that the postulated linearity was convincingly re- 
covered from the ML analysis performed. 

In the B and V colors the slope of the relationship is 
the same while in R and in /, in particular, it significantly 
differs. In the literature the inverse value of a is usually 
given. Using the value of a obtained for B would give a bias 
of about 0.44 log[iVo(m)/iV(m)] mag in the estimation of 
the optical extinction, in the / color. 

3.4. Surface distribution of the obscuring material 

Based on the calibration procedure one may assign extinc- 
tion to each pixel in the star count maps, in all of the four 
colors studied. Fig. [3 shows the contour maps of the ex- 
tinction obtained in this way. The main body of the cloud 
stands out with a contrast of several magnitudes from the 
less obscured region behind the bow shock. 

The overall distribution of the obscuring material, ob - 
tained from our study and that of iKandori et alJ l|2003l) . 
has a reasonably good correlation. There are, however, re- 
markable differences between them, in particular in the 
densest part of the cloud. The probable reason for these 
discrepancies lies in the different ways used in obtaining 
extinction maps from the surface distribution of the stars 
in the region of LI 251. Both studies had a mesh of 2' reso- 
lution but applied different kind of smoothing. We used a 
boxcar of 6'x 6' size while the other study was a Gaussian 
filtering of 6' resolution. It is probably important to note 
that the tails of the Gaussian filter give contributions from 
much wider areas, in particular in the densest part where 
the star counts have low values. 

A further reason for the discrepancy might originate 
from the method of converting the star count maps into 
extinction. Our method accounts for the foreground stars 
whose contribution increases the counts and decreases cor- 
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5 



FIELD 


AREA{SqD) 


as 


eras 


av 








ai 




1. 


0.9980 


0.45 


0.05 


0.30 


0.05 


0.12 


0.05 


0.10 


0.03 


2. 


0.8770 


1.55 


0.10 


1.25 


0.07 


0.85 


0.07 


0.63 


0.07 


3. 


0.1515 


3.30 


0.15 


2.65 


0.15 


2.20 


0.15 


1.75 


0.2 


4. 


0.0605 


6.5 


0.80 


5.2 


0.80 


> 5.3 




4.5 


0.4 



Table 3. Summary of the star count analysis of L1251, distance moduli; the weighted mean of the data gives 7.58 ±0.2 



FIELD AREAjSqD) tb (Trg ry a,-^, m cr,.^ n 

1. 0.9880 7^5 060 >8.50 '- '- ~ '- ~ 

2. 0.8770 7.50 0.35 7.75 0.55 <7.6 - <7.4 

3. 0.1515 7.75 0.30 7.65 0.45 7.65 0.45 <7.4 

4. 0.0605 7.05 0.7 <7.4 - <7.4 - <7.3 




Fig. 2. Contour maps of the extinction in the B, V, R and / color. The 'flying bullet' form of the main body and the 
bow shock displayed in Fig. ^is clearly visible. The coordinates of the (0,0) position are same as in Fig. 1. 
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respondingly the estimated extinction v alues. We made 
simulations of the lWainscoat et aD l)l992|) model assuming 
a limiting magnitude of about 19 mag and 6 mag of extinc- 
tion at the distance of L1251. The conventional method 
of star count extinction conversion resulted only in a 5 
mag. value, i.e. one magnitude less than the true (6 mag.) 
extinction. 

4. Discussion 

4.1. Shape of the cloud 

We have already indicated above that the extinction map 
derived from the star counts resembles a body flying at 
hypersonic speed across an ambient medium. 

Th e shape of LI 251 was accounted by iToth et "aD 
lll995^ for a shock wave passing the cloud and produced 
by a nearby supernova. They showed that the cooling by 
the H2 molecules plays an important role in the formation 
of the observable shape of the cloud. 

The presence of a bow shock as indicated by 
the extinction maps suggests another type of cloud- 
environment interactions. The blunt body solutions of hy- 
personic flows are standard topics in textbooks (see e.g. 
[Haves fc Probstein (,19591 . Fo r the astrophysical context 
of the probl em see the works of [l^pzvczka fc Tenorio-Tasfl3 
(TgS?) and ' Canto & Ragal l|l998l) l. One can identifv the 
tail of L1251 with a wake of the head of the cloud, typical 
of blunt body hypersonic streaming patterns. 

In the following we do not try to fit the form of the 
bow shock using a solution of the blunt-body problem. 
This solution would allow a better understanding of the 
role of different significant physical parameters defining 
a particular fit, however this is beyond the scope of the 
present paper. An obvious significant parameter would 
be the Mach number of the flow which is easy to cal- 
culate from the angle between the tw o asymptotes of the 
bow shock I Haves fc Probsteiiilll959fl . The Mach number 
is given by M = sin(a/2)~^ where a is the angle between 
the asymptotes. It yielded M « 2 in our case. 

The symmetry axis of the bow shock has a tilt of about 
10 deg to the main axis of L1251. The symmetry axis of 
the bow sh ock points toward th e center of the bub ble dis- 
covered bv lGrenier et all l|l989l) . iKim et all l(2n0(t found 
a runaway star which also might have been caused by the 
SN explosion that was probably responsible for the bubble 
and derived an age of 10^ yrs. According to this picture 
the bow shock resulted from an encounter of the cloud 
with the wind coming from the interior of the bubble. 



4.2. Properties of the obscuring material 

The estimation of the total extinction in different colors 
makes it possible to get the value of the total to selective 
extinction Ry = av/Es-v ■ The scatterplot between ay 
and as is displayed in Fig.|3| The relationship between ay 
andas can be written in the form of as = (l-|-l/i?y)xay. 
We marked with lines (labelled with the corresponding 



Rv values) in Fig. 3 the relationships between ay and 
ttB assuming that Ry = 3 (close to the canonical value 
for the general interstellar matter) and Ry — 6. The line 
assuming that the interstellar absorption does not depend 
on the color (grey approximation) is also marked with 
Ry = inf. 

Figure 13 clearly demonstrates that above ay = 3 (i.e. 
within the main central body of LI 251) the points depart 
from the Ry — 3 line and approach the lines of higher 
Ry as one moves to higher extinction. Departure from 
the canonical Ry = 3.1 towards higher values is a com- 
mon b ehavior of dense clouds (see the review bv iMath'i^ 
l|l99(t V 

Postulating a relationship in the form of = ax 
ay, where a = (1 -I- 1/Ry), and substituting Ry values 
between 3-6, we can infer that a depends only weakly 
on ay and may have approximately a form of a{ay) = 
ao+ai xay. Assuming that Ry = 1 + 1/aQ ~ 3.1 we made 
a least squares fitting to get ai. In this way we obtained 
Ry{ay) = 1/(0.32 - 0.019 x ay) and as = a{ay) x ay 
gave a reasonably good fit to the points in Fig. |21 

Based on this result we may assign a value of Ry to 
each pixel in the ay extinction map of L1251. We obtained 
a map of the total to selective extinction as given in Fig. 
One may infer from this map that the high values of 
Ry are concentrated only in the densest parts of L1251. 

The canonical dependence of the interstellar extinction 
on the wavelength is re alized by the standard relation of 
ISavas-e Ma.thi j l|l 979h which can be well represented by 
a linear function of the extinction on 1/A in the range of 
the B, V, R, I colors used in our study. 

Departure from the standard Ry — 3.1 value above 
ay = 3 in dense interstellar cloud s, as we found in the 
case of L1251, was also obtained bv lWhittet et al liool) 
in the Taurus region. They claimed that for extinctions 
ay > 3 real changes in grain properties may occur, char- 
acterized by observed Ry values in the range of 3.5-4.0. 
A simple model for the development of Ry with ay sug- 
gested that Ry may approach values of 4.5 or more in the 
densest regions of the cloud. According to IWhittet et alJ 
(2001) the transition between "normal" and "dense cloud" 
extinction occurs at ay « 3.2, a value coincident with the 
threshold extinction above which H20-ice is detected on 
grains within the cloud. 

The Ry values derived in our analysis correspond to 
those obtained by iKandori et al.l l|2003l) within the lim- 
its of the statistical errors. There are differences, however, 
in the surface maps between the two studies. To obtain 
Ry values we used the relationship between the as and 
ay values so our surface map reflec ts that of the e xtinc- 
tion in the V color. On the contrarv. lKandori et al.l l(2003f) 
applied an adaptive averaging of B — V and V — I color 
indices by varying the size of the smoothing window, keep- 
ing constant the number of stars in it but at the cost of 
the spatial resolution, in particular in the densest part of 
the cloud. 
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L1 251 : surface distribution of R R 




X offset (arcmin) 



Fig. 3. Scatterplot between ay and (left panel). The relationship between ay and as can be written in the form 
of Ob — (1 + 1/^) X fly- We marked with lines (labelled with the corresponding R values) the relationships between 
av and as assuming that R ^ 3 (close to the canonical value for the general interstellar matter) and R — 6. The 
line assuming that the interstellar extinction does not depend on the color (grey approximation) is also marked with 
R = inf. Assigning Ry to each pixel we obtained a map of the total to selective extinction (right panel). One may 
infer from this map that the high values of Ry are concentrated only in the densest parts of L1251. The coordinates 
of the (0,0) position are same as in Fig. 1. 



Table 5. Mass of the cloud (using the formula of lDicknianl 

(dzi) 



FIELD 


AREA (SQD) 


MASS (Mq) 


NO. OF FIX. 


1. 


0.9880 


153.34 


887 


2. 


0.8770 


490.48 


784 


3. 


0.1515 


212.60 


136 


4. 


0.0601 


132.05 


54 


5. 


0.0089 


26.07 


8 



4.3. Mass of the cloud 

Knowing the distance of L1251 we converted the ay ex- 
tinction values into the mass using the empirical formula 
given by Dickman (1978). According to this formula 

M = {adf^l—Yay{i) (8) 
ay ^ 

where M, a, d, /i, ay(i) are the mass of the cloud, the an- 
gular size of a pixel, the distance, the mean molecular 
mass and the extinction of a pixel, respectively; Nr j ay — 
1.87 X W'^cmr'^mag-^ and Nh = Nhi + 2Nh2 - Based on 
Dickman formula we computed the mass of the different 
subregions of the cloud as given in Table [S] 

Fields 3, 4 and 5 represent the main body of L1251. 
Based on Table |S1 the total mass of this part of the cloud 
is 371 Mq. Thi s value c an be compared with the figure 
obtained by iSato et alJ l)l994() based on C^^O measure- 
ments, taking into account that the field occupied by our 
main body lies completely inside the region covered by the 
C^*0 study and contains only 85 % of the mass calculated 



from it. Re-scaling this fraction of mass to the distance we 
obtained in this paper one gets 422M0, surprisingly close 
to our estimated value. Studying a somewhat larger area 
than ours, Lec (1994) obtained 610 from his ^^CO and 
^^CO measurements. 

4.4. Mass model of the head of L1251 

In the following we try to model the head of the cloud 
assuming spherical symmetry. 

Using the da ta obtained with VLT of the ESO 
lAlves et al.l ll200lh modelled the spatial structure of the 
Bok globule Barnard 68. They suppo sed the globule to be 
an is othermal Bonnor-Ebert sphere l)Bonnoilll956t lEbertl 
19551) in pressure equilibrium with the external, much hot- 
ter and less dense interstellar medium. 

We have not required isothermality in our analysis. 
The isothermal solution is a specific case of the class o f 
polytropic spheres studied extensively by lEmdenl lll907l) . 
The isothermal solution represents the case when the poly- 
tropic index n = oo. Actually, the mass of an isother- 
mal sphere is infinite and the limit between the finite 
and infinite mass solutions is n = 5. Although the mass 
is finite in this case the size of the sphere is infinite. 
One obt ains finite size solutions only in the n < 5 case. 
Recently, iMedvedev fc Rvbickil l|200l|) studied the proper- 
ties of polytropic spheres near n = 5 and found them very 
suitable to characterize the structure of molecular cloud 
cores. 

We fixed the mass center of the sphere at the maxi- 
mum value of ay. This choice is quite reasonable in view 
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Fig. 4. Fitting the radial distribution of the surface mass 
density by a Schuster sphere. Up to 2.5 pc distance from 
the center the fit is very good. Beyond this distance, how- 
ever, there is an excess of the observed mass density due 
to the tail of the cloud. 

of Equation (jSJ . We then averaged the ay values over con- 
centric annuli. This procedure gave a radial profile of the 
extinction of the head of LI 251. This profile can be con- 
verted into the surface mass density by Equation (|5]). 

Evaluating the radial profile obtained from the ob- 
served data we projected the n = 5 polytropic sphere 
llSchustedllSsl to get its surface mass density distribu- 
tion. The Schuster sphere has two free parameters to be 
adjusted: the central mass density and a scale parameter. 
We get the fit displayed in Fig. 01 Up to 2.5pc distance 
from the center the fit is very good. Beyond this distance, 
however, there is an excess of the observed mass density 
due to the tail of the cloud which significantly distorts the 
spherical symmetry of the head. 

The density parameter of the fitted Schuster sphere, 
4:5Mqpc~^ (3.06 X lO^^^gcm^'^), gives an estimate of the 
central mass density of the head of L1251. Comparing the 
finite size polytropic solutions of n < 5 with those of the 
Schuster sphere (i.e. n — 5) one can infer that they con- 
centrate more mass at finite distances and give a much 
worse fit to the points observed. 

Assuming a polytropic gas sphere one can compute 
the radial velocity dispersion profile from the fitted den- 
sity profile, based on the polytropic equation of state. The 
p/p ratio of the pressure and mass density gives the T 
temperature, i.e. the velocity dispersion. The velocity dis- 
persion projected onto the celestial sphere can be directly 
compared with that observed. The projected FWHM of 
the V-profile of a Schuster sphere is displayed in Fig. 
El along with tho s e of t he NH^ molecule measured by 
Eoth fc WalmslevI lIlQQfit) . 

There is a considerable scatter of the measured line 
widths in the head of the cloud around the mean which 
is well matched by the Schuster solution. The Schuster 
curve, however, gives an unexpectedly good fit in the tail 
region. This good fit is not expected far from the head 
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Fig. 5. FWHM of the projected velocity profile of a 
Schuster sphere (dashed hne) and the NH^ molecule mea- 
sured bv lToth k WalmslevI l)l996|) . There is a considerable 
scatter of the line widths in the head of the cloud around 
the mean which is well matched by the Schuster solution. 
The Schuster curve gives an unexpectedly good fit in the 
tail region. 

since due to the tail the mass distribution drastically dif- 
fers from the spherical symmetry. However, we may as- 
sume that our predicted FWHMs of th e NHS l ine refe rs to 
the in itial stage and those measured bv lToth'fc WalmslevI 
|)199(tI) to the final stage of the formation of density en- 
hancements in the tail region. We conclude therefore that 
the formation of structures in the tail left the linewidth 
practically unchanged. This indicates an isothermal con- 
traction. However, the question remains open: what kind 
of instability played a significant role in forming the den- 
sity enhancements in the tail? We note that according to 
iToth fc WalmslevI lll99 6h the thermal energy is dominant 
over the turbulent energy in the NHS cores of the tail re- 
gion and thus we cannot exclude the scenario that thermal 
instability played a significant role in the cloud fragmen- 
tation. 

5. Summary and conclusions 

In this paper we studied the _B, V, R, I star count maps 
of L1251 to derive the spatial structure of the obscuring 
matter in the cloud. We assumed that the pixel values in 
the star count maps form a one dimensional manifold in 
the BVRI four dimensional parameter space and defined 
areas of similar extinction by means of multivariate k — 
means clustering. After defining areas of equal extinction 
we derived the amount of obscuration using a maximum 
likelih ood procedure based on a Monte Carlo simulation 
of the IWainscoat et al.l l)l992|) model. As a byproduct of 
the model fitting we obtained the SSO ± SOpc distance of 
the cloud. 

The extinction maps clearly showed the main body of 
the cloud and a less dense region having a form of a bow 
shock of a blunt body. The form of the bow shock allowed 
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us to calculate the approximate Mach number {M « 2) of 
the streaming around the head of L1251. 

Comparing the obscuration of the dust in different col- 
ors we calculated the dependence of the total to selective 
extinction {Ry value) on the visual extinction. The total 
to selective extinction exceeded Ry = 5 in the densest 
part of the cloud. 

Using Equation (|SJ to convert the optical extinction 
distribution into surface mass density we obtained the 
mass of 371 Mq for the cloud, in reas onably good agree- 
ment with the 410 Mq obtained bv ISato et all (|B94), 
based on C^^O measurements. 

Assuming a spherical symmetry for the head of L1251 
we computed the radial distribution of the mass and com- 
pared it with a polytropic model of n = 5 (Schuster 
sphere). Up to r = 2.5 pc from the center of the mass 
the fit is excellent but beyond this distance the observed 
points start to depart remarkably from the fit due to the 
drastic distortion of the spherical symmetry by the tail of 
L1251. 

In the head of LI 251 the Schuster model fit matches 
well the mean FWHM of the 1.3 cm NH3 line measured by 
iToth fc WalmslevI l)l996l) . It gives an unexpectedly good 
fit in the tail region indicating that isothermal contraction 
played a significant role in forming the density enhance- 
ments. 

Acknowledgements. This work was partly supported by OTKA 
grants T 024027, T 34584 and T 037508 of the Hungarian 
Scientific Research Fund. P.A. acknowledges the support of 
the Bolyai Fellowship. 



— . 1998, ApJS, 115, 59 

Kun, M. & Prusti, T. 1993, A&A, 272, 235 

Kun, M., Vinko, J., & Szabados, L. 2000, MNRAS, 319, 777 

Lee, C. W., Myers, P. C, & Tafalla, M. 1999, ApJ, 526, 788 

— . 2001, ApJS, 136, 703 

Lee, Y. 1994, Journal of Korean Astronomical Society, 27, 159 
Mathis, J. S. 1990, ARA&A, 28, 37 
Medvedev, M. V. & Rybicki, G. 2001, ApJ, 555, 863 
Murtagh, F. & Heck, A. 1987, Multivariate data analysis 

(Astrophysics and Space Science Library, Dordrecht: Reidel, 

1987) 

Nikolic, S., Johansson, L. E. B., & Harju, J. 2003, A&A, 409, 
941 

Rozyczka, M. & Tenorio-Tagle, G. 1985, A&A, 147, 220 
Sato, F. & Fukui, Y. 1989, ApJ, 343, 773 

Sato, F., Mizuno, A., Nagahama, T., et al. 1994, ApJ, 435, 279 
Savage, B. D. & Mathis, J. S. 1979, ARA&A, 17, 73 
Schuster, A. 1883, British Assoc. Report, 427 
Sume, A., Downes, D., & Wilson, T. L. 1975, A&A, 39, 435 
Toth, L. v., Horvath, A., & Haikala, L. A. 1995, Ap&SS, 233, 
175 

Toth, L. V. & Walmsley, C. M. 1996, A&A, 311, 981 
Wainscoat, R. J., Cohen, M., Volk, K., Walker, H. J., & 

Schwartz, D. E. 1992, ApJS, 83, 111 
Whittet, D. C. B., Gerakines, P. A., Hough, J. H., & Shenoy, 

S. S. 2001, ApJ, 547, 872 
Zhou, S., Wu, Y., Evans, N. J., Fuller, G. A., & Myers, P. C. 

1989, ApJ, 346, 168 



References 

Alves, J. ., Lada, C. J., & Lada, E. A. 2001, Nature, 409, 159 
Balazs, L. G., Eisloeffel, J., HoU, A., Kelemen, J., & Kun, M. 

1992, A&A, 255, 281 

Benson, P. J. & Myers, P. C. 1989, ApJS, 71, 89 
Bonnor, W. B. 1956, MNRAS, 116, 351 
Cambresy, L. 1999, A&A, 345, 965 
Canto, J. & Raga, A. 1998, MNRAS, 297, 383 
CaseUi, P., Benson, P. J., Myers, P. C, & Tafalla, M. 2002, 
ApJ, 572, 238 

Cernicharo, J., Bachiller, R., & Duvert, G. 1985, A&A, 149, 
273 

Cudaback, D. D. & Heiles, C. 1969, ApJ, 155, L21 

Dickman, R. L. 1978, AJ, 83, 363 

Dieter, N. H. 1973, ApJ, 183, 449 

Ebert, R. 1955, Zeitschrift Astrophysics, 37, 217 

Emden, V. R. 1907, Gaskugeln (Leipzig) 

Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, P. C. 

1993, ApJ, 406, 528 

Grenier, I. A., Lebrun, F., Arnaud, M., Dame, T. M., & 

Thaddeus, P. 1989, ApJ, 347, 231 
Hayes, W. D. & Probstein, R. F. 1959, Hypersonic Flow Theory 

(Academic Press New York and London) 
Kandori, R., Dobashi, K., Uehara, H., Sato, F., & Yanagisawa, 

K. 2003, AJ, 126, 1888 
Kendall, M. & Stuart, A. 1976-1979, The Advanced Theory of 

Statistics, Vols. 1-3 (Macmillan) 
Kun, M. 1982, Astrophysics, 18, 37 



